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SUMMARY 


A diamond wing and body combination was designed to have an area 
distribution which would result in near optimum zero-lift wave-drag 
coefficients at a Mach number of 1.00, and decreasing wave-drag coeffi- 
cients with increasing Mach number up to near sonic leading-edge condi- 
tions for the wing. The airfoil sections were computed by varying their 
shapes along with the body radii (blending process) to match the selected 
area distribution and the given plan form. The exposed wing sections had 
an average maximum thickness of about 3 percent of the local chords, and 
the maximum thickness of the center-line chord was 5*^9 percent. The wing 
had an aspect ratio of 2 and a leading-edge sweep of 45°. Test data were 
obtained throughout the Mach number range from 0.20 to 3 • 50 at Reynolds 
numbers based on the mean aerodynamic chord of roughly 6,000,000 to 
9 , 000 , 000 . 

The zero-lift wave-drag coefficients of the diamond model satisfied 
the design objectives and were equal to the low values for the Mach number 
1.00 equivalent body up to the limit of the transonic tests. From the peak 
drag coefficient near M = 1.00 there was a gradual decrease in wave-drag 
coefficient up to M = 1.20. Above sonic leading-edge conditions of the 
wing there was a rise in the wave -drag coefficient which was attributed in 
part to the body contouring as well as to the wing geometry. 

The diamond model had good lift characteristics, in spite of the 
prediction from low-aspect -ratio theory that the rear half of the diamond 
wing would carry little lift. The experimental lift-curve slopes obtained 
at supersonic speeds were equal to or greater than the values predicted by 
linear theory. Similarly the other basic aerodynamic parameters, aero- 
dynamic center position, and maximum lift -drag ratios were satisfactorily 
predicted at supersonic speeds . 
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INTRODUCTION 


The design of airplanes is largely a matter of compromise. This is 
indicated by theoretical considerations of the selection of a wing plan 
form as reported in references 1, 2 , and 3* For instance, the diamond 
wing which is the best structural shape of the zero-taper-ratio plan form 
is generally considered aerodynamic ally the poorest. Aerodynamic is ts 
point out on the basis of low-aspect -ratio or slender-body theories, that 
the rear half of a low-aspect -ratio diamond wing does not produce lift 
(ref. 1 or 4) , and the zero sweep of the mid-chord line results in high 
transonic wave drag (refs. 1 through 5)« For most low-aspect -ratio wing- 
body combinations, the wave drag at Mach number 1.00 can be reduced a 
large amount by the application of the simple area rule of reference 6; 
however, at Mach numbers only slightly above 1.00 the wave-drag coeffi- 
cients will generally increase considerably (in particular for config- 
urations with unswept wings). Similarly the supersonic area rule of 
reference 7 permits one to design for supersonic Mach numbers, but the 
procedures are more complicated and will result in an increase in wave- 
drag coefficients at nondesign Mach numbers. The area rules are usually 
applied by selecting a wing of fixed plan firm and thickness distribution 
and shaping a fuselage to provide a smooth area distribution. A more 
flexible approach used in the present investigation was to vary both the 
wing thickness distribution and the fuselage shape to obtain a blended 
wing -body combination. 

The purpose of the investigation was to seek the lowest possible wave 
drag at Mach number 1.00 (configuration mad^ optimum at M = 1.00), which 
was consistent with a decreasing supersonic wave -drag coefficient. A 
brief statement of the design concept involved is to concentrate the wing 
volume near the body center line with a favorable or minimum distortion 
with Mach number of the area distributions lefined by the theory of refer- 
ence 7* Details of the concepts and procedures used in the design of the 
test model (aspect -ratio -2 diamond wing) are given in appendix B. The 
symbols used in the report (including appendix B) are listed in appen- 
dix A« 

Experiments were conducted on the diamond wing -body combination at 
Mach numbers from 0.20 through 3*50 and at Reynolds numbers per foot that 
ranged from 1,000,000 to 6,000,000. Data were taken at each Mach number 
at a Reynolds number per foot of at least 3*000,000 to ^-,000,000. The 
Reynolds numbers based on the mean aerodynamic chord of 26.67 inches would 
be, of course, over twice the numbers stated. Also tested at transonic 
speeds was a body with the same area distribution as the wing-body combi- 
nation. As a check on the experimental results, theoretical predictions 
were made of the friction drag, zero-lift wave drag, drag due to lift, 
lift -curve slope, maximum lift-drag ratio, and aerodynamic center position 
of the wing -body combination. 
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MODELS AND TESTS 


The geometric details of the diamond wing and body combination are 
presented in figures 1 through 3 and in tables I and II. The radii for 
the body components are listed in table I* The forward and rearward body 
components are defined by the design area distribution, KA 25 , shown in 
figure 2. The central portion of the body was designed to be as small as 
practical with sufficient space for the strain-gage balance which measured 
the aerodynamic forces* The low-drag shape selected for this central 
portion was part of a von Karman ogive (l - 4o in* and r-^ - 1.625 in.) 
defined by the following equation (also see table i): 
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and was a cylinder for the body stations rearward of 40 inches. ^ Experi- 
mental wave -drag results and pressure distributions for a von Karman 
ogive are reported in reference 8. The M = 1.00 equivalent body of the 
wing -body combination was also tested and is defined in table I and 
figure l(b). 

The wing coordinates are listed in table II. The wing -thickness 
distribution was computed as described in appendix B and is illustrated 
in figure 3- The wing thickness in each case is formed by straight-line 
elements perpendicular to the model center line forming triangular span 
vise sections. Representative maximum thicknesses in percent chord are 
5 . 49 , 4.04, and 3.07 for center -line, body juncture (2 in. from the 
center line), and mean aerodynamic chords respectively. The average 
maximum thickness of the exposed wing was about the same as the maximum 
thickness of the mean aerodynamic chord. 

Both the diamond model and its equivalent body were tested in the 
Ames l4-foot transonic wind tunnel at Mach numbers of 0.60 through 1.20. 
The diamond model was also tested in the Ames 12-foot subsonic pressurized 
wind tunnel and in the 9 - by 7 -foot and the 8 - by 7 -foot supersonic test 
sections of the Ames Unitary Plan wind tunnel. The latter supersonic 
facilities are described in reference 9 * A photograph of the diamond 
model in the 8 - by 7 -foot supersonic test section is shown in figure 4. 

The equivalent body is shown in the perforated test section of the 
transonic tunnel in figure 5 . The Mach number range of the tests in each 
facility as well as the range of the other test variables is shown in the 
following table: 
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Wind tunnel throat 

M 

R/ft 

a, 

deg 

Trans it ion 

12-foot 

0.20 

3,000,000 aid 

-4 to 26 

Free 

12 -foot 

0.50 

6,000,000 

<3 

0 
-p 

1 

Free 

l 4 -foot 

o.6o to 0.80 

3,500,000 to 

-2 to 11 

Free and fixed 



4 , 000,000 



l 4 -foot 

0.80 to 1.20 

0 

8 

*\ 

0 

0 

0 

«\ 

-V 

-2 to 6 

Free and fixed 

9- by 7“foot 

1-55 

4,000,000 

-2 to 7 

Fixed 

9“ by 7“foot 

1.55 to 2.35 

3,000,000 

-2 to 11 

Fixed 

8“ by 7 “foot 

2.50 to 3*00 

3,000,000 and 

-2 to 12 

Free and fixed 



4,000,000 



8- by 7 “foot 

3.00 to 3.50 

3,000,000 

-2 to 12 

Free and fixed 

8- by 7~foot | 

2.50 to 3-50 

1,000,000 

-2 to 12 

Free 


Three -component aerodynamic forces and moments were measured and 
corrected by standard procedures. For the model sizes and shapes, the 
force corrections for blockage and buoyancy were generally found to be 
negligible. Wall interference corrections wore required for the angle 
of attack and drag data obtained in the subsonic wind tunnel, and these 
corrections were made based on the theory of reference 10 . At all Mach 
numbers, the drag coefficients were adjusted by equating the body base 
pressures to free-stream static pressures. All aerodynamic coefficients 
are based on the complete plan-form area of the diamond wing of 800 square 
inches. The pitching -moment coefficients were computed about a longitu- 
dinal center 34*50 inches rearward of the nose of the body or 7*5 percent 
of the root chord forward of the centroid of the wing area. This position 
was selected for neutral longitudinal stability at supersonic speeds. 

The procedures used to compute the friction drag and to fix transi- 
tion, and the effects of the boundary-layer :rip are discussed in 
appendix C; however, a brief description wil_ be given here of the dis- 
tributed roughness used to fix the location of the bound ary -layer transi- 
tion from laminar to turbulent flow. The grit size required (mean height 
of about 0.040 in.) was determined by prelim. nary tests at M = 3-50 and 
was used on the body and the wing at all Mach numbers although the size 
was excessive for the lower Mach numbers. A laminar-flow area of 5 percent 
of the wing area was allowed (i.e., the grit was located 1.13 in. rearward 
of the wing leading edge and of the body nose in a streamwise direction). 




RESULTS AM) DISCUSSION 


The basic data are presented in figures 6 through 10. These figures 
illustrate that the effects of the Reynolds number changes were generally 
small; however, the Reynolds numbers will be noted where the effect on 
the drag data may be of some importance. The discussion will deal mainly 
with data taken at a Reynolds number per foot of 3,000,000 or 4,000,000. 
These basic data plots also indicate that the effect of fixing transition 
is relatively small and primarily effects the drag data. The effect of 
fixing transition is discussed in appendix C, so the following discussion 
is directed toward transition-fixed results. 

The discussion is presented in two parts: the first section is an 
analysis of the zero-lift wave drag (the major design parameter), and 
the second is an analysis of the aerodynamic trends with Mach number. 


Zero -Lift Wave Drag 


A fundamental way of experimentally evaluating how well the diamond 
model met the design objective (to seek a low wave drag at Mach number 
1.00, which was consistent with a decreasing supersonic wave-drag coeffi- 
cient) is to compare its wave drag with that of its equivalent body which 
is in turn compared with an optimum body. It should be remembered that 
the wave drag of a wing-body combination is usually greater than that for 
■Its M = 1.00 equivalent body at Mach numbers only slightly above 1.00. 
Only transonic data were obtained for the equivalent body and thus the 
comparisons are made at these speeds. At transonic speeds the diamond 
model fully met the design expectations as shown by the wave-drag coeffi- 
cients in figure 11. Essentially identical zero-lift wave-drag coeffi- 
cients were obtained for the diamond model and Its equivalent body at 
Mach numbers up to 1.20. From the peak drag coefficient near M - 1.00 
there was a gradual decrease in drag coefficient up to M = 1.20. 

In figure 12 the wave-drag coefficients of the equivalent body are 
compared with those for a minimum wave -drag body and with theoretical 
predictions based on the method of reference 11. The data points are 
shown for direct evaluation. The coefficients are based on the plan-form 
area of the diamond wing. Each body had a closed-body fineness ratio of 
12-5 and the body shapes are discussed fully in appendix B. The data, 
for the Sears -Haack body (minimum transonic wave-drag body for prescribed 
volume and length, ref. 12) were obtained from the investigation of 
reference 13- Both bodies had transition fixed; however, the grit used 
on the Sears -Haack body was selected for transonic speeds and was much 
smaller (about 0.003 in. in diameter). Note that the wave-drag coeffi- 
cients for the equivalent body are substantiated by theory and are as 
low as or lower than the minimum wave-drag body for Mach numbers greater 
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than 1.04. The difference in wave -drag coefficients at M = 1.00 between 
the bodies is much greater than that indicate! by theory (see appendix b) 
and may be partially due to the difference in the size of the grit used 
to fix transition. 

The diamond model was designed for M = 1.00 and one would not 
expect the wave -drag coefficients to continue to decrease through Mach 
numbers at which the velocities normal to the wing leading edge are sonic 
or supersonic. It is of interest to see how experiment and predictions 
would compare at higher supersonic Mach numbers. Such comparisons are 
made in figure 13 ) with theoretical predictions based on 4-0 harmonics 
representing the derivatives of the area curves. Predictions were also 
made with 25 harmonics; however, the larger number of terms was required 
to give a good representation of the area distributions. At transonic 
speeds the solutions for the diamond model based on 25 harmonics were 
satisfactory and almost equal to the 49-term solution. 

The theoretical computations of wave drag were made by two methods, 
each of which was based on the procedures of reference 11. Each method 
would give essentially the same answer if the body alone equivalent area 
distribution did not vary with Mach number (i.e., slender bodies accept- 
able to the theory (ref. 7))- Thus at the higher supersonic Mach numbers 
where different answers were obtained as shown in figure 13 , both theo- 
retical predictions are questionable. It is of interest to note that 
this fact plus an analysis of the area distributions indicate that the 
rise in wave -drag coefficients at Mach numbers near sonic leading-edge 
conditions of the wing may be attributed in pjirt to the body contouring. 
Method 1, used for the results presented in figure 13, is the more 
convenient of the two methods, and it involve; the use of various combi- 
nations of equivalent area distributions obtained from oblique cuts of 
the total wing-body combination as a function of roll angle 0 . It is 
theoretically more correct to handle the variation of the body area 
independently of the wing area, because for a body of revolution the area 
distribution varies only with Mach number and not with roll angle 0 . 

This latter procedure was used for the result; listed in figure 13 as 
method 2 . Experimental data were not available at sonic -leading -edge 
conditions; however, it is evident that both theoretical methods tend to 
overestimate the wave drag under these circumstances. Note in figure 13 
that at a Mach number of 3 . 00 the higher Reynolds number data agree better 
with the theory. The results of method 1 agree better with the experi- 
mental results at the higher supersonic speeds , so these theoretical 
results were used in determining the theoretical lift -drag ratios presented 
in the next section of the report. 
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Aerodynamic Trends With Mach Number 


Presented in figure l4 are summary plots of experimental and theo- 
retical values of drag due to lift, aerodynamic -center location, lift- 
curve slope, and maximum lift-drag ratios for the diamond model. At 
supersonic speeds the lift-curve slope and the aerodynamic -center position 
were computed with equations of linearized theory available in refer- 
ence 1 or l4. The theoretical results at subsonic speeds were computed 
from the curves of reference 15 which is based on a simplified lifting- 
surface theory. The theory at M = 1.00 is based on low-aspect -ratio 
theory (ref. 1 or and is directly applicable only to delta wings. 

Except for the zero-lift drag data, all the theoretical results are for 
the wings alone without thickness. The friction drag was computed as 
outlined in appendix D and the wave drag at zero lift was computed by 
method 1 as discussed previously. The maximum lift-drag ratios were 
predicted from a combination of these results. 


The theoretical drag due to lift parameter 



as shown in 


figure l4(a), was computed as the reciprocal of the theoretical lift- 
curve slope at all Mach numbers and therefore was without leading-edge 
thrust. Also shown in figure l4(a) is the value, = l/itA, for 

minimum vortex drag due to lift. The results indicate that the sharp- 
edged diamond model did have some leading-edge thrust at subsonic and 
transonic Mach numbers; however, for Mach numbers greater than 2.50 the 
experimental drag due to lift parameters approached the theory without 
leading-edge thrust, and theory and experiment were in good agreement. 

The subsonic tests were without the transition fixed, but the laminar- 
flow region for the low-speed tests should be comparable to the 5 percent 
allowed (1.13 in. streamwise). 


The variation in aerodynamic -center position with Mach number 
occurred primarily at transonic speeds as shown in figure l^-(b) . The 
subsonic and supersonic levels were reasonably well predicted, although 
the theoretical shift in aerodynamic -center position was less than that 
obtained experimentally. 


Probably the most interesting result of this investigation is the 
good lift characteristics of the diamond plan form, in spite of the 
prediction that the rear half of the wing would carry little lift accord- 
ing to low-aspect -ratio theory. At M = 1.00 the experimental lift- 
curve slope (fig* l^-(d)) for the diamond model (0.0596) was higher than 
the theoretical value for a similar aspect-ratio-2 delta wing (0.05^8). 
Similarly, for Mach numbers equal to 1*55 or greater (where data were 
available) the experimental lift -curve slopes were equal to or greater 
than the values predicted by linear theory (see fig. l^(d)). 
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The maximum lift -drag ratios of figure l4(c) reflect the previously 
discussed zero-lift wave -drag coefficients, the drag due to lift, and the 
lift-curve slopes. The excellent transonic zero-lift wave-drag coeffi- 
cients, some leading-edge thrust, and high ] if t -curve slopes all helped 
to produce high maximum lift -drag ratios of 9*5 or better at transonic 
speeds . At the higher supersonic speeds the agreement between theory and 
experiment is very good, but probably would be improved even further for 
a thinner wing and a body with less contouring. 


CONCLUDING REMARKS 


The zero -lift wave -drag coefficients of the diamond model were 
essentially identical with the low values for its M = 1.00 equivalent 
body up to the limit of the transonic tests (M = 1.20). From the peak 
wave-drag coefficient near M = 1.00 there was a gradual decrease in wave- 
drag coefficient up to M = 1.20. Near sonic leading-edge conditions of 
the wing there was a rise in the wave-drag coefficient which was attributed 
in part to the body contouring but was primarily attributed to the wing 
geometry. 


The diamond model had good lift characteristics, in spite of the 
possibility that the rear half of the diamond wing would carry little 
lift according to low -aspect -ratio theory. The experimental lift -curve 
slopes obtained at supersonic speeds were equal to or greater than the 
values predicted by linear theory. Similarly the aerodynamic center posi- 
tion, and the maximum lift-drag ratios were satisfactorily predicted at 
s upe r s on ic s pe e ds . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., July 8, 19*9 
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APPENDIX A 
SYMBOLS 


A 

% 

An 

b 

2 

C D 




aspect ratio 

area distribution of the approximate Newtonian nose shape 

coefficients determining the magnitude of the harmonics of a 
Fourier sine series 

semispan 

drag coefficient 

(All aerodynamic coefficients are based on the total wing 
area. ) 

slope of the curve of drag coefficient due to lift vs. lift 
coefficient squared., taken at the lift coefficient data point 


nearest to that for 




PA 


max 


zero-lift drag coefficient 



C R 

c 

E» (m) 
k 



lift coefficient 
lift-curve slope, per deg 

pitching-moment coefficient about body station 34.50 measured 
from the body nose 

center line, or wing-root chord 

mean aerodynamic chord 

complete elliptic integral of the second kind 
average height of transition grit 

M = 1.00 equivalent body radii or area distribution 
maximum lift -drag ratio 


l 


model length 



Mach number 


p cot A 

total number of terms used in computing A C-n 

specific term or harmonic used to define an area distribution 

Reynolds number 

body radius 

body maximum radius 

body base radius 

cross-sectional area 

total wing area 

wing semithickness 

body station, measured aft from model nose 

aerodynamic center location, where x x = station measured from 
leading edge of wing-root chord 

volume 

volume relative to a Sears -Haack body with minimum wave drag for 
given volume and length (ref. 12) 

angle of attack 

n/m 2 - 1 

zero -lift wave-drag coefficient 

roll angle of a cutting plane tangent to a Mach cone as measured 
between the z axis and the intersection of the cutting plane 
with the yz plane 

angle of sweepback of leading edge of wing 

cotangent of sweep angle of plan-form leading edge 
cotangent of sweep angle of plan-form trailing edge 

transformation of the length x to radians, cos” 1 


l 
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APPENDIX B 
DESIGN PROCEDURE 


As mentioned in the introduction to this report, the general approach 
was to seek the lowest possible wave drag at Mach number 1.00 which was 
consistent with a decreasing supersonic wave-drag coefficient. The basic 
design concept was to concentrate the wing volume near the body center 
line with a favorable or miniinuin distortion of a selected area distri- 
bution with changes in Mach number. The distortions in the area distri- 
butions are the result of obtaining equivalent area distributions from 
oblique cuts tangent to the Mach cone according to the theory of 
reference 7* 

The diamond plan form selected was not considered as ideal from a 
wave-drag standpoint, but rather as a probable shape of interest to air- 
plane designers. However, by selecting a plan form with zero taper ratio, 
the wing volume would naturally tend to be concentrated near the body 
center line. The symmetry of the diamond plan form justified its selec- 
tion, because the distortions of its area curves could also be made 
symmetrical. 

With the plan form specified the next step was to select a spanwise 
thickness distribution. The thickness distribution selected was formed 
by straight-line elements forming triangular spanwise sections (prepen- 
dicular to the model center line), which incidentl y results in a rigid 
structural shape. Various cusp shapes were also considered which would 
further concentrate the wing volume along the model center line or improve 
the supersonic wave drag. Along this line of reasoning an idea of R. T. 
Jones is simply illustrated in sketch (a), where it is shown that the 



Sketch (o) 
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spanwise thickness distribution can be altered for some higher design 
Mach number by adding zero area to the M = 1.00 area distribution. This 
idea was not investigated very fully,, but an application with parabolic 
arcs, rather than the straight lines shown in sketch (a), was tried for 
the diamond plan form for Mach numbers up to 1.20. Little improvement in 
computed wave drag was indicated over the straight-line thickness distri- 
bution; so these more complex distributions were not justified for the 
conditions considered. For higher Mach numbers or for other configu- 
rations, complex spanwise thickness distributions might be justified. 

Other interesting types of wing thickness distributions formed by straight 
lines are discussed in reference 1 6 . 

The area distribution was next selected to help compensate for the 
distortion of the wing area distribution with Mach number. The area 
distribution was also selected to have an equivalent body with near opti- 
mum theoretical zero-lift wave drag. The area curves considered in the 
selection of the design area distribution are shown in figure 15- The 
length and the maximum cross-sectional area were held constant, although 
a reduction in volume was allowed relative to the Sears-Haack body (mini- 
mum transonic wave-drag body for given volume and length, ref. 12). The 
area-distribution curve for this Sears-Haack body is shown in the lower 
left-hand corner of figure 15, and the volume ratio is of course equal to 
1.00. The derivative of this area curve is completely defined by the 
amplitude of the second harmonic of a Fourier sine series. It may be 
noted that all the other bodies had a reduction in volume and generally 
an increase in drag. A reduction in zero-lift wave -drag coefficient was 
computed for one area curve, defined by two terms of a sine series, shown 
in the upper right-hand corner of figure 15- The upper left-hand curve 
is that for a three-quarter power body, an approximation of the Newtonian 
nose shape (based on impact theory, ref. 17 ) > reflected to form a closed 
body. The theoretical zero-lift wave-drag coefficients are quite similar 
considering that the one area curve requires almost an infinite number of 

terms (N = 25 Tor Z^C-n, shown) for definition while the other curves 

^o 

required only one or two terms. The coefficients are based on a wing 
area of 800 square inches. 

A preliminary analysis of these area distributions with the diamond 
plan form was made to indicate how the theoretical wave-drag coefficients 
would vary with p cos 6 or roughly Mach number as follows : 
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Plan 

Area distribution 
(see fig. 15 ) 

Drag coefficient, ^D 0 “ 25 ) 

form 


A 3 

a 4 

\l >4 

p cos 0=0 

p COS 0=0.6633 

P cos 0=3.3541 

Diamond 

% 

a n 

a w 

a n 

0.0034 

0.0023 

o.oo 4 o 

Diamond 

1.696 

0 

0 

0 

.0025 

.0059 


Diamond 

1.405 

0 

-.727 

0 

.0026 

.0053 


Diamond 

1-551 

0 

-.360 

0 

.0023 

.0053 



■''Sears -Haack body, minimum wave drag for given volume and length from 
linearized theory. See reference 12 . 


The p cos 0 value of O.6633 represents the maximum cutting angle for 
M = 1.20 and the value of p cos 0=0 represents all cutting angles for 
M = 1 . 00 . From these data it was indicated that the Aj\j area distri- 
bution with the diamond wing fully met the desired design objectives, 
because the wave-drag coefficient is near optimum at M = 1.00 and 
decreases with increased Mach number up to 1 . 20 . For the highest value 
of p cos 0 (3.35^1) the theory was not relied upon for the preliminary 
analysis, since the velocity components normal to the wing leading edge 
were supersonic, however, the ind.icated wave-drag coefficient was reason- 
ably low. Thus it was decided to use the (Ajj) area distribution; however, 
some modifications to the selected area distribution were made, as shown 
in figure l6. To be consistent with the drag coefficients shown in 
table V, the Ajj distribution was replaced by the distribution defined 
by 25 harmonics (A 25 ). Next the maximum cross-sectional area was corrected 
to the original value (to maintain a fixed maximum cross-sectional area) 
as shown in figure l6 to result in the KA 25 area distribution which 
was the design area curve used. This area distribution was cut off at 
station 60, with a curve with zero slope at the base and defined by 25 
harmonics, to permit sting -mounting of the models. These modifications 
to the design area curve resulted in only a small change in the computed 
wave-drag coefficient from the original value of 0.0034 to O.OO366 and 
O.OO369 for N = 25 and N = 49 , respectively. 

With the area distribution and the wing spanwise thickness distri- 
bution specified (triangular spanwise sections) the wing section could be 
computed directly if the model was a wing alone configuration. For the 
present investigation, forward and rearward body components were required 
because the selected length of the area curve was greater than the wing 
streamwise length. Also a minimum body size was specified in the central 
region of the model to permit installation of the strain-gage balance. 

For an actual airplane, similar body components would be a function of 
the engine size or cockpit size. Temporarily ignoring the minimum body, 
the next step in the model design was to compute the thickness distri- 
bution of the wing along the model center line as shown in sketch(b). 
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Bod y stati on x 
Sketch (b) 


The wing thickness at the leading and trailing edges of the root chord 
approached infinity as the span approached zero. In the region where the 
minimum body radii were greater than half the wing thickness, the combined 
area distribution would be greater than that specified and thus the exposed 
minimum-body area was removed from the given area distribution and the 
wing thickness was recomputed. The convergence was rapid because of the 
small size of the exposed minimum-body area. The blending process was the 
last step in the design. The root section of the wing was made to have 
zero thickness at the leading and trailing edges, and the thickness of 
the wing and the radii of the body were varied or blended as shown in 
sketch (c) to match the given area curve and to avoid abrupt changes in 



surface slope in the wing or the body components. The body blending of 
course did not occur unless the wing extended outside the body. The 
center line or root airfoil sections had the greater surface slopes and 
the outboard sections were almost biconvex as may be seen in figure 1 and 
interpreted from figure 3- The blended regions of the final body are 
also noted in table I. 
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APPENDIX C 


FIXING BOUNDARY -LAYER TRANSITION AND FRICTION DRAG 


Confidence in the results of a zero -lift wave -drag investigation 
can be obtained only if the variations in friction-drag coefficients with 
Mach number are known or may be computed accurately. This is because the 
experimental wave drag is determined by subtracting the estimated friction 
drag from the measured zero-lift drag. The laminar and turbulent friction 
drag were computed at M = 0 with zero heat transfer by the method of 
reference 18. The variation of the friction-drag coefficients with Mach 
number was computed from the following equations : 

Laminar flow, reference 19 

C Df ^ 

(c Df ^) M =o 

Turbulent flow (smooth surface), reference 20 


Turbulent flow (rough surface), reference 21 


= (l + r M 2 ^j (3) 


where the recovery factor, r = 0.86 and 7 = 1.1-. This latter equation 
is required for the wing area covered by a distributed roughness. 

Because the extent of the laminar boundary -layer flow will vary with 
Reynolds number and Mach number, a rigorous wave -drag comparison between 
models can be obtained only if identical regions of laminar, turbulent, 
and turbulent -rough flow are maintained at each test condition. Fixing 
transition at high Mach numbers and low Reynolds numbers on wings is 







difficult as indicated in reference 22. The single -element roughness 
used in reference 22 resulted in large drag (considered to he primarily 
wave drag) and was not adequate in fixing transition at M = 3*53 and 
Reynolds numbers per foot of approximately 1,400,000 to 2,700,000. This 
appendix described what was done to fix transition on the diamond model 
at M = 3-50 and R per foot = 3*000,000, aid what effect fixing transi- 
tion had on the aerodynamic data. 


SELECTION OF DISTRIBUTED ROUGHNESS 


Prior experience rejected single -element roughness for the boundary- 
layer trip and suggested the use of a distributed roughness. Based on 
estimates of the short length of the laminar -flow at subsonic speeds at 
R per foot = 4,000,000 the following grit location and extent was selected. 
The distributed roughness was formed by first spraying a thin coat of 
cement and then spraying (at low pressure) the grit on a strip of the wing 
0.22 inch wide and 1.13 inches rearward of the wing leading edge and the 
body nose (measured in a streamwise direction). This resulted in fixing 
the laminar region to 5 percent of the wing area and the rough region to 
1 percent of the wing area. Grit concentration was approximately 200 per 
square inch and this value was used in the computation of friction drag. 

The size of the grit was determined from preliminary tests in the 3- by 
Y-foot supersonic wind tunnel. The end of the laminar-flow region was 
determined by the use of fluorine sublimation material as discussed in 
reference 23* The grit finally selected had a mean height of 0.040 inch 
and varied from about 0.030 to 0.050 inch. The grit location on the wing 
and the grit size is indicated in the full-scale photograph of figure 17- 
In figure 17 the painted white lines are located 3 inches apart in the 
streamwise direction. Laminar flow is indicated ahead of the grit, mixed 
vortex and laminar flow are indicated by tie parallel streaks, and fully 
turbulent flow is indicated where the exparding black streaks coalesce. 

In this preliminary test the grit was not placed close enough together 
in a couple of regions as indicated by the extended white streaks repre- 
senting localized regions of laminar flow. More uniform grit spacing 
was used when the force data were obtained* 

An example where transition was not considered to be fixed is shown 
in figure l8(a) for a Reynolds number per loot of 1,000,000 (the white 
fluorine crystals were not removed from the wing tip). The satisfactory 
performance of the 0.040-inch grit at R per foot = 3,000,000 and M = 3 . 50 
is shown for comparison in figure 18 (b)* I he grit was sufficiently large 
that it extended out of the boundary layer; however, if the grit were 
located further rearward, the length of the laminar-flow region would be 
increased at high supersonic Mach numbers and would be decreased at 
subsonic Mach numbers. Smaller sized grit in the same location, of a 
size to lie within the boundary layer, did not fix transition. These 
results indicate a need for further investigation of the parameters used 
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in the suggested roughness Reynolds numbers recommended in references 24, 
25, and 26. In particular the use of the curves of reference 2 6 for 
Mach numbers of the order of 3. 00 or greater does not apparently apply 
for cases such as the present investigation where transition is desired 
near the wing leading edge. The roughness Reynolds number for the case 
of figure l8(b) is about 7,000 or 5 to 7 times greater than that indicated 
by the curves of reference 2 6 for two-dimensional and. three-dimensional 
flow, respectively. Perhaps the velocity and temperature used in the 
roughness Reynolds number should be based on other boundary -layer param- 
eters such as the momentum or displacement thickness which would better 
group the results for roughness particles both within and partially 
outside the boundary layer. The sublimation procedure used in the present 
investigation to indicate transition (with the criteria that there should 
be no laminar flow of consequence rearward of the grit) is probably more 
conservative than the detection of spots of turbulence used, in refei 
ences 24, 25 , and 26. 


EFFECT ON THE AERODYNAMIC DATA 


Previously presented data of figure 7 indicated that fixing transi- 
tion had very little effect on the lift or moment data and the greatest 
effect on the zero-lift drag coefficients. A summary plot of the zero- 
lift drag coefficients for the diamond model is presented in figure 19* 

The addition of the distributed roughness increased the zero-lift 
drag coefficients at the higher supersonic Mach numbers to a lesser degree 
than at subsonic Mach numbers (ACp 0 = 0.0013 to 0.0017> respectively). 

Also the experimental increase in zero-lift drag is very nearly the same 
as the computed increase in friction drag (note, in particular, the higher 
Reynolds number data at M = 3.00). These two facts qualitatively indicate 
that the grit did not produce wave drag in spite of its large size. The 
increase in drag due to fixing transition with the grit was a reasonable 
amount and was estimated quite well as shown in figure 19. With transi- 
tion free tests, the line of transition between laminar and turbulent flow 
as indicated by the sublimation material was quite irregular at subsonic 
speeds so the uncertainty in predicting the laminar region is indicated 
by the shaded region. The base -drag coefficients which were typical of 
all bodies tested with or without wings are included in this figure for 
completeness of all drag components involved at zero lift, and illustrate 
that the body base may contribute to the zero-lift drag if the base is 
not completely filled with jet exhaust. 

Since the effect of fixing transition was not very large on the zero- 
lift drag coefficients, the effect on the maximum lift-drag ratios was 
not large as shown in figure 20. The greatest reduction in (L/p) may 
was about 0.8 at subsonic speeds, and only 0.2 to 0.4 at supersonic speeds. 
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TABLE I.- COORDINATES FOR BODIES, INCHES 


M = 1*00 

equivalent body 

ka 25 

X 

r 

0 

0 

.375 

.075 

.750 

.148 

1.500 

.258 

2.250 

.356 

3.000 

.449 

3.750 

.542 

4.500 

.623 

5.250 

.693 

6.000 

.761 

6.750 

.830 

7.500 

.903 

8.250 

.976 

9.000 

1.048 

9-750 

1.116 

10.500 

1.179 

n.250 

1.238 

12.000 

1.296 

12.750 

1.353 

13.125 

1.388 

13.875 

1.443 

14.625 

1.509 

15-375 

1.570 

16.125 

1.630 

16.875 

1.688 

17.500 

1.733 

18.000 

1.762 

19.000 

1.823 

19-875 

1.896 

20.625 

1.948 

21-375 

2.000 

22.125 

2.054 

23-625 

2.166 

24.375 

2.221 

25-125 

2.275 

26.625 

2.374 

27.375 

2.420 

28.125 

2.464 

29.625 

1 2.555 

30.375 

31.125 

2.603 

2.655 


Diamond 

model 


KA pi5 body 


19-360 

20.000 

21.000 

22.000 

24.000 

24.800 

26.000 
26.400 
28.000 
29.600 
31.200 

32.800 
33.6OO 


Blended 

body 


1.860 
1.874 
1.825 
1-712 
1.471 
1.392 
1 . 3^7 

Von Karman I. 36 O 
ogive I- 405 

1.447 
1.487 
1.522 
1-539 


M = 1.00 
quivalent body 
KA 25 


32.625 
33-375 

34.125 

35.625 

36.375 
37.500 

38.625 

39-375 

40.125 

41.625 

42.375 

43.125 

44.625 

45.375 

46.125 

46.875 

47.625 

48.000 

48.375 

49.125 

49.875 

50.625 
51-375 

52.125 

52.875 

53.625 

54.375 
55-125 

55.875 

56.625 
57.300 
57 - 600 
57.900 
58.200 
58.500 
58.800 
59.100 
59 - 400 
59.7OO 
60.000 


2.765 
,820 
2.873 
2.957 
2.984 
3.000 
2.984 
2-957 
2.919 
2.820 
2.765 
2-709 
2.603 
555 
_ 509 
2.464 
2.420 
2.4-00 
2.374 
326 
2.275 
2.221 
l66 
2.110 
2.054 
,000 
,948 
I.896 
1.846 
1.795 
1.747 
1.728 
1.709 
1.692 
1.677 
1.662 
1.649 
1.637 
1-629 
1.625 


Diamond 

model 


34.400 Von Karman l-55 lf 


3^.200 ogive 

36.000 
36.8OO 
37 . 600 
38.400 
39.200 

40.000 


569 
1.582 
1.594 
1.605 

1.6l4 

1.621 

1.625 


Cylinder 


51.000 

52.000 

53.000 

54.000 

55.000 
55.640 


Blended 

body 


1.625 

1.657 

1.730 

1.827 

1.878 

1.860 


KA a5 body 
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TABLE II.- COORDINATES FOR DIAMOND WING, INCHES 





Figure !•- Sketch of the model. 
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A-23952 

Figure 5*~ Equivalent body (KA 25 area distribution) in the test section of the Ames l4-foot 

transonic wind tunnel* 










Tronsition fixed 
Transition free 



Figure 7.- Transonic aerodynamic characteristics of the diamond model with transition fixed and 
free as determined from tests in the Ames l4-foot transonic wind tunnel. 
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Figure 7 • “ Continued* 
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Figure 7.- Concluded. 





Figure 8.- Supersonic aerodynamic characteristics of the diamond model with transition fixed as 
determined from tests in the 9~ Uy 7"foot supersonic test section of the Ames Unitary Plan 
wind tunnel# 
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Figure 8.- Continued. 





Figure 8.- Concluded 




Figure $>•“ Supersonic aerodynamic characteristics of the diamond model with transition fixed and 
free as determined from tests in the 8- by 7“foot supersonic test section of the Ames Unitary 
Plan wind tunnel. 
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Figure 9*" Continued 
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Figure 9 *“ Continued 





Figure 9 *“ Continued 


















Figure l4.- Aerodynamic trends with Mach number for the diamond model as indicated by experiment 

and theory (transition fixed except as noted). 












miodli _ sinin + IJ^i , /2x 

n-l ntl J ’ cosv \ l 1 



sponding theoretical (N = 25) zero-lift wave-drag coefficients are for M = 1.00. 
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Figure l6*- Derivation of the design area distribution 
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A-23660 

Figure 17*- Full-scale photograph of the 0.040 inch grit used to fix transition on the diamond 
wing (preliminary test at M = 3.0 and Ji/ft = 3,000,000). 
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(a) M = 3.5, R/ft = 1,000,000, k = 0.060 inch 



A-23675. 1 

(h) M = 3 . 5 , R/ft = 3,000,000, k = 0.040 inch 


Figure l 8 .- Conditions for which the fluorene sublimation material indi 
cated unsatisfactory and satisfactory fixing of the boundary -layer 
transition. 









